Electron Spin Dynamics in Impure Quantum Wells for Arbitrary Spin-Orbit Coupling 
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Strong interest has arisen recently on low-dimensional systems with strong spin-orbit interaction 
due to their enhanced efficiency in some spintronic applications. Here, the time evolution of the 
electron spin polarization of a disordered two-dimensional electron gas is calculated exactly within 
the Boltzmann formalism for arbitrary couplings to a Rashba spin-orbit field. It is shown that the 
classical Dyakonov-Perel mechanism of spin relaxation gets deeply modified for sufficiently strong 
Rashba fields, in which case new regimes of spin decay are identified. These results suggest that 
spin manipulation can be greatly improved in strong spin-orbit interaction materials. 
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PACS numbers: 72.25.Rb, 72.25.Dc 71.70.Ej 

The physics of transport of electron spins in low- 
dimensional systems has become a popular theme of re- 
search due to the possible impact in future electronic 
applications 0, 0. Key subjects of studies concern the 
problem of controlling the electron spin polarization 
through the tailoring of the spin-orbit (SO) interaction 
0, 0] , and the knowledge of the physical parameters gov- 
erning the spin relaxation time t s . Main sources of SO 
coupling are the Rashba interaction arising from struc- 
tural inversion asymmetries of low-dimensional struc- 
tures such as quantum wells or wires |5(, and the Dres- 
selhaus interactionpresent in bulk crystals lacking sym- 
metry of inversion^. In the presence of Rashba and/or 
Dresselhaus interactions, t s basically arises from the ran- 
domization of the SO (pseudo)magnetic field induced by 
momentum scattering with impurities or other interac- 
tions (Dyakonov-Perel mechanism 0). 

Of great interest for spintronic applications are mate- 
rials with strong SO interaction since they are more ef- 
fective spin-manipulators Q and spin-current generators 
through the spin-Hall effect 0. In this respect, systems 
of interest may be, for example, HgTe quantum wells 
whose Rashba interaction leads to SO band splitting A so 
as large as ~ 30 meV0, or surface states of metals and 
semimetals which display giant SO band splittings of the 
order of 0.2 eV or surface states of metals and semimet- 
als which display giant SO band splittings of the order 
of 0.2 eV or larger [ToL Ull. or even the heavy Fermion 
superconductor CePt 3 Si [13 whose lack of bulk inver- 
sion symmetry leads to A so ~ 50 — 200 meV 0. For 
all these systems, the SO splitting A so is no longer the 
smallest energy scale, as in most semiconductors, and can 
be comparable to the Fermi energy or even larger. 

Despite that the time evolution of the spin polariza- 
tion S has been thoroughl y st udied in the weak SO limit 
A so /e_F = H EH El, onl y few partial studies 
have been devoted to the strong SO case A so /ep ^ 
[rH Hil . Of particular interest is the problem of as- 
sessing the role of A so /e^ on the spin relaxation mech- 
anism and of clarifying to which extent the classical 
Dyakonov-Perel (DP) description gets modified by fi- 
nite values of A so /ep. 



In this letter, a generalized kinetic equation for S is 
formulated for arbitrary strengths of the SO interaction 
and in the presence of spin-conserving coupling with im- 
purities. For the case of a two-dimensional electron gas 
confined in the x — y plane subjected to the Rashba in- 
teraction, the kinetic equation for the z-component of S 
is solved explicitly for any value of A so /ep. It is found 
that for a sufficiently strong Rashba interaction the DP 
relaxation mechanism gets modified, with the spin po- 
larization displaying a slow (power law) decay with time. 
Furthermore, in the extreme A so /eF 3> 1 limit, a fast ex- 
ponential decay is obtained with t s proportional to the 
momentum scattering time r p , i.e., spin relaxation is en- 
hanced by disorder. These findings suggest that tailoring 
of spin memory through A so / 1 tp can be much more effec- 
tive than previously thought. 

Let us consider an electron gas whose non-interacting 
hamiltonian Hq is: 



H 



„C ks 



(1) 



where c fc (c ks ) is the creation (annihilation) operator for 
an electron with momentum k and spin index s =1,1, a 
is the spin- vector operator with components (a x ,a v ,a z ) 
given by the Pauli matrices, and Ok is a k dependent SO 
pseudopotential vector whose explicit form is not essen- 
tial for the moment. In the above expression, e k is the 
electron dispersion in the absence of SO coupling. Let us 
consider as momentum-relaxation mechanism a coupling 
Vkk' to spin-conserving impurities located at random po- 
sitions R.;. Additional interaction channels as those pro- 
vided by phonons or electron-electron couplings can be 
treated as well by following the same steps described be- 
low. Assuming that Vkk' is switched on at time t = to, 
the time evolution of S(i) for t > to can be obtained from 
the equation of motion of the density matrix p(t) 20]: 

= -j[EM'{t) - pkk'W^k'] - ^fkk'W, (2) 
at n n 

where /Skk' (t) is a 2 x 2 matrix with elements jO ks ,k's' (i) = 



2 



(ks|/>(t)|kV), £? k = e k - fj, + fflk-CT, and 

fW(i) = £ [v kp e iRi -( k -P)p pk ,(t) 



-l/ pk ,e lR -(P- k VkpW 



(3) 



Equation © can be formally integrated and, after the 
usual adiabatic (t — > — oo) and Markov approximations 
HHEl Pkk'(0 can be expressed as: 

/>00 

(t)~ / dt'e- l5t 'e^ l - gkt '/' i r kk ,me^' t ' /?i (4) 



Pkk 



where the limit S — > + must be taken after the inte- 
gration. By using the anticommutation property of the 



Pauli matrices, a{5j + UjCJi = l<5y , the exponential op- 
erators in Eq.© can be put in a form more suitable for 
integration over t': 



z±l E k t'/n = 1^2 e ±lE ^'/ h (l + aCl w ■ a), (5) 



where a — ±1, fj k = O k /|f2 k |, and E k ± = e k ± ||f2 k | 
are the eigenvalues of Hq [Eq.©]. At this point, averag- 
ing over the impurity positions Ri restores the transla- 
tional invariance: (p kk ' (i))i, np = <5 k ,k'PkW; an d> finally, 
the equation of motion of electron spins for a given wave 
vector k, S k = |Tr(er/9 k ), is obtained by integrating over 
t' . By retaining only the scattering contributions (Boltz- 
mann approximation) |l9| . the final result is therefore: 



dt 



tt k x s k - Yl V X\ E {t 1 ~ a ^ ■ ^p)]( s k - s P ) + a$n k [n p ■ (s k - s p )] + apn^ ■ (s k - s p )] 

P a/3 

+h( a n k + pn p )(fr - / p )} 5{E ka - E p/3 ). (6) 



r 



where n is the impurity concentration and / k = ^Tr(p k ) 
is the electron occupation distribution function. Equa- 
tion © is very general, and holds true for both bulk 
and low-dimensional systems with no restrictions on the 
particular form of J~i k . In this letter, Eq.© is used 
to find the time evolution of the electron spin polar- 
ization S = ^ k S k for a two-dimensional electron gas 
confined in the x — y plane and subjected to a Rashba 
SO interaction. In such a case, the effective SO field is 
S^k = {inky, —Jnk x ,0) and the corresponding SO split 
bands depend solely on the magnitude k of the wave vec- 
tor and can be written as E-^ a = E^a = ^—{k + cekpY' \ 
where m is the electron mass and fc^ = -^Ir |22j . 
In the following, instead of using A so = hy^kp, the 

SO splitting will be parametrized by the Rashba energy 
h 2 k 2 

€ rt = 2m = it^R' that is the minimum inter-band exci- 
tation energy for an electron sitting at the bottom of the 
lower band [see Fig©a)] . 

By assuming that the momentum dependence of the 
impurity potential can be neglected, Vk p = V, S k is con- 
veniently decomposed in its even, S k , and odd, S k , parts 
with respect to k. It is then clear that S = J^k S k and, 



from Eq. © , dS/dt = £ k ^k x S k 



J2k ^k x S k , where 



f2 k = Jl k has been used. Of course, dS/dt is also 
equal to ^2kdS^_/dt, so that, after a further derivative 
with respect to time, the equation of motion can be re- 
cast in the following form: 



E 



dt 2 



o dS k 

dt 



0. 



(7) 



By taking the vector product with Jl k , the odd part of 
Eq.© can be rewritten as: 



dS£ 

dt 



(tt k xS k )- 


2imV 2 
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^ka — Ep/s), 





a/3 p 



(8) 



where the summation over momenta has cancelled all 
terms odd in p and the identity (f2 k ■ ri p )(fi k x S°) — 

(ft k x^ p )(fi k -s p ) - [n p xs p -n k -(si p xs°)n k ](k/p) 

has been used. From Eq.©, the equation of motion of 
S k is instead given by: 



dt 



O k xS k 



2irnV 2 
Ah 



EE(* 

a/3 p 



k B p 



-ficrfi k /k - h/3n p f°) 8{E ka - E. 



(9) 



where /£ = |(/ k - /_ k ) is the odd part of / k . 

By using Eqs.© the term O k x dS^/dt in Eq.© can 
be eliminated in favor of S k and S k . Next, by using 
Eq.©, also the terms containing S k can be eliminated 
and Eq. © reduces to an equation of motion of the com- 
ponent S k only, that is sufficient to find the time evolu- 
tion of S = J2k S k - Let us consider the z-component of 
S, S z , for which the presence of / k in Eq.© has not effect 
since these terms have zero component in the z-direction. 



3 



In this way, Eq.Q reduces to: 



°°dkk 
2tt 



d 2 S\ 



zk 



dt 2 



2 ; 2 



Xk \ q C 
^1) ^ 



T k dS 



zk 



dt 



= 0, 



(10) 



where S% k = f Q * d(f> S'| k /27r, with 4> being the angle be- 
tween the directions of k and the x-axis, r" 1 = ^-nV 2 N 
is the momentum relaxation rate for a two-dimensional 
electron gas with zero SO splitting and density of states 
(DOS) iV = m/27rfi 2 , and: 



a/3 



^rJ2l ^P k -T p )5{E ka ~E pf3 ). (12) 

A solution to Eq. (|10|) is obtained by equating to zero the 
expression within square brackets, which leads to a ho- 
mogeneous differential equation of the second order for 
S zk (t). In this way the functions T k and Xk assume re- 
spectively the meaning of renormalization of the damp- 
ing term and of shift of the (bare) precessional frequency 
7_r/c. It can be easily realized from Eas. (|llll2fl that in the 
weak SO limit ef/cr —> oo, for which E k ± — > both 
the damping renormalization and the frequency shift are 
absent (1^ = 1 and Xk = 0), indicating that these quanti- 
ties stem from additional intra- and inter-band scattering 
channels opened when ep/cR is finite. Let us take a closer 
look at T k and Xfe by performing the integration over p 

Hum 



dpp 



in Eqs. 



k R — k 



0<k<k 



R 



T k = 



1 + P k *- k \ k * k R <k< 2k R 



{k—kn.)k 



(13) 



Xk 



2k 2 +3k R k 

k 2 — k 2 
(2fe«-fc)fc 

(k-2k R ) 2 



I o 



(k-k R )(3k R -k) 

(4k R -k)k R 
(k-k R )(k-3k R ) 



2k R < k 

< k < k R 
k R <k< 2k R 
2k R <k < 3k R 
3k R <k < Ak R 
4k R < k 



(14) 



For k > 4fcfl, T k and Xk are the same as in the zero SO 
limit, while for lower momenta they acquire a strong k 
dependence (plotted in Fig^c)) arising from the com- 
bined effect of the reduced dimensionality (D — 2) and 
the SO interaction. In particular, the divergence of T k 
at k = k R and those of Xk at k = k R and k — 3k R 
are due to scattering processes probing the SO induced 
van Hove singularity of the DOS of the lower sub-band 
which diverges as N^(E) cx d e R /E [see FigQIb)]. As it 
is shown below, such strong k dependence has important 
consequences on the spin polarization dynamics. 

Let us now turn to evaluate the explicit time depen- 
dence of S z . From Eq. (fTUfl , a general solution for 5| fc (t) is 




0k1 2 k 1 2 3 4 1 
' k/k R * N ± (E)/N 



2 3 
k/k 



FIG. 1: (a): Rashba SO split electron dispersions ^-(fcztfcn) 2 
in units of 6r — ^-"(r- The horizontal dashed line indicates 
the Fermi level, (b): density of states for the ± bands, (c): 
plots of r fe , Eq.l^J, and Xh, EcLltTIl. 



t ) , where 



given by a linear combination of exp^- — 

S fc = T' 2 . — Xk — (^ T plRk) 2 , whose coefficients are fixed 
by imposing some initial conditions. If at t = elec- 
trons have been prepared with a non-equilibrium spin- 
state occupation but equilibrium distribution for each 
spin branch then, at the lowest order in the initial weak 
spin imbalance Sfi — (//j — /i|)/2, 



S e zk (0) is simply: 



^ e fc (o) 



Sfi 

lR.k 



^2 af (E ka - /i), 



(15) 



where fo(x) — (exp(x/T) + 1) _1 is the Fermi dis- 
tribution function, T is the temperature and p, > 
is the chemical potential. Furthermore, by imposing 
that limt^oo IS^^)! < oo and by arbitrarily choosing 
dS z (0)/dt — for r" 1 = 0, at zero temperature (/i = e^) 
S z (t) is readily found to be: 



S,(t) =- 



5)i 

2-K^R 



dk 



0(E fc -r^)exp - 



2r n 



+0(T 2 k - S fc )exp( -t-M coshf 



2r n 



2r n 



(16) 



where 8 is the unit step function, fc^ = k R {\J £f/cr ± 1) 
for 6f/€r > 1 and k^ = k R (l ± \/7fJer) for e F /e R < 1 
[see Fig^a)]. For cf/cr ™ > oo, fc^ — > kp and Eq. i|TT)|) 
reduces to the classical formula 



(17) 

where S* z (0) = —HSfiNQ and Q R = j R kp is the Rashba 
frequency which characterizes the (damped) spin pre- 
cession behavior S z (t) ~ S z (0) exp(— ^-) cos(il R t) for 
2r p Q, R ^S> 1 and the DP relaxational decay S z (t) ps 
S z (0)exp{-T p n R t) for 2t p Q r < 1 0. 

For finite values of ep/e R , Eq.^J) starts to deviate 
from the classical regime (|17|) . Let us first consider 
£f/cr > 25. In this case, the integration over k in 
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FIG. 2: (a): Spin polarization time evolution from Eo. 11611 
(solid line) for t p = 3.3 fs, en — 5 meV and ef = 25 meV. 
The dotted line is the power law decay of Eg. l20H while the 
grey solid line is Eg. 11711 . (b): Crossover from the power law 
decay of 1201 and the fast relaxation decay of (I2H obtained 
for r„ = 3.3 fs, £f = 5 meV and several values of tR. 



Ea. (|16fl spans values necessarily larger than 4fc#, where, 
according to Eas. i|13ll4l) . = 1 and Xk — [see also 
FigJIJc)]. For weak impurity scattering (2r P 7^A:< ^> 1) 
Efe is negative and Eg. (j!6|) reduces to: 



S z (t) «S*(0)exp 



t \ sin(f2>t) — sin(fi<t) 



2r T 



(o> - o<)t 



(18) 



where 



-fii k^. The main feature of Ea. fl8p is that 
S z (t) oscillates with two different Rashba frequencies 
associated with the two SO splitted Fermi surfaces, 
FigUJa). Two distinct Rashba frequencies characterize 
also the relaxation regime obtained when the scattering 
with impurities is strong enough that 2r p 7flfc> 1 (and 
so Efc > 0) holds true. Also in this case the integration 
of Eq. i|lb| is elementary and 



^jrerf^TpfiS.i) - erf (y^ 

~2 y/Tp^-Q^t 



(19) 



where erf is the error function. 

The spin precession and relaxation regimes of 
Eas. l|18ll9|) are governed solely by the enhanced mo- 
mentum phase space settled by finite values of €f/^r- 
Instead, for cf/cr < 25, also the momentum depen- 
dence of Tfe and %fe becomes relevant, leading to im- 
portant anomalous features of the spin dynamics. One 
of these is particularly striking and it is found when 
1 < ep/cR < 9. In this case, the integration over k in 
Eq. (|Lo|) crosses the point k = 2kn where Xk changes sign 
[see Fig^c)]. Hence, if 4:T p jnkji -C 1, V^fc can be ex- 
panded as Tfc — (2kn — k)/3kpi for k < 2kR, while when 
k > 2kn \k becomes negative leading to exponentially 
small contributions to S z (t) for sufficiently long times. 
Eq. (|16fl then can be approximated to: 



S z (t) 



5fi 



2 k 



]R dk (i*R-*) t 
— e 6r p fc i? 
2 



35 z (0)r p 



(20) 



The surprising result of Ea. H20(l provides the rather in- 
teresting prediction that, for sufficiently strong SO inter- 
action and momentum scattering, the spin polarization 
decays as a power law rather than exponentially. In this 
case therefore the memory of the initial spin polarization 
can be much longer lived than in the DP regime, as shown 
in Fig. Ufa). Another striking feature is that obtained 
in the extreme cf/cr <C 1 limit in which the integration 
of Ea. lfTojl becomes restricted to a narrow region around 
k = kn where both Tk and Xk diverge as l/\k — Hr\, so 
that Ea. (frp|) becomes: 



S z (t) oc exp 



8r c 



(21) 



indicating that for extremely strong SO interaction, mo- 
mentum scattering increases the spin polarization decay. 
The power decay of l|20|l and the fast relaxation regime of 
(|21|l are plotted in Fig|2Jb) from a numerical integration 
of Eq.|(TSJ| for r p = 3.3 fs, ep = 5 meV and €r ranging 
from 50 meV down to 1 meV. 

To conclude, the kinetic equations describing the time 
evolution of the spin polarization have been formulated 
for arbitrary strength of the SO interaction. Explicit 
solutions for quantum wells with Rashba-like SO interac- 
tions predict the failure of the DP relaxation formula for 
sufficiently strong SO couplings. In particular, the mem- 
ory of the initial spin state can be strongly enhanced or 
reduced depending on €f/cr, suggesting an alternative 
route for spin manipulation in spintronic applications. 

I thank E. Cappelluti for fruitful discussions. 
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